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ACTIVE GALACTIC NUCLEI SELECTED FROM GALEX SPECTROSCOPY: THE IONIZING SOURCE 
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ABSTRACT 

We use a complete sample of Lya emission-line selected active galactic nuclei (AGNs) obtained 
from nine deep blank fields observed with the grism spectrographs on the Galaxy Evolution Explorer 
( GALEX) satellite to measure the normalization and spectral shape of the AGN contribution to the 
ionizing background (rest- frame wavelengths 700-900 A) at 2 ^ 1. Our sample consists of 139 sources 
selected in the redshift range z — 0.65 — 1.25 in the near-ultraviolet (NUV; 2371 A central wavelength) 
channel. The area covered is 8.2 deg^ to a NUV magnitude of 20.5 (AB) and 0.92 deg^ at the faintest 
magnitude limit of 21.8. The GALEX AGN luminosity function agrees well with those obtained using 
optical and X-ray AGN samples, and the measured redshift evolution of the ionizing volume emissivity 
is similar to that previously obtained by measuring the GALEX far-ultraviolet (FUV; 1528 A central 
wavelength) magnitudes of an X-ray selected sample. For the first time we are able to construct the 
shape of the ionizing background at 2 ~ 1 in a fully self-consistent way. 

Subject headings: cosmology: observations — galaxies: active — galaxies: distances and redshifts — 
intergalactic medium 



1. INTRODUCTION 

An important parameter in the cosmological model- 
ing of galaxies and the intergalactic gas is the ionizing 
background radiation produced by the overall popula- 
tion of active galactic nuclei (AGNs) and star-forming 
galaxies as a function of wavelength and redshift. Usu- 
ally the input ionizing source spectrum for such models is 
computed by convolving an approximate quasar spectral 
energy distribution (SED) determined from composite 
quasar spectra with a quasar luminosity function deter- 
mined from optical data (e.g., Haardt & Madau 1996; 
Madau et al. 1999; Haardt & Madau 2001; Meiksin 2005; 
Siana et al. 2008; Dall'aglio et al. 2009). (Note that 
Gilmore et al. 2009 adopted the Hopkins et al. 2007 
quasar luminosity function, which tries to take into ac- 
count obscured sources, but obscured sources will not 
contribute to the ionizing background.) However, using 
optical quasar luminosity functions that have been ex- 
trapolated beyond where they have been measured (es- 
pecially if the turn-over has not been observed) can re- 
sult in overestimates of the ionizing background. X-ray 
luminosity functions, which probe to fainter luminosities 
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(e.g., Barger et al. 2005; Richards et al. 2005), provide a 
better approach (Cowie et al. 2009; hereafter, CBT09). 

In CBT09 we measured the ionizing fluxes from a wide- 
field (0.9 deg^) X-ray sample with optical and UV imag- 
ing observations and found that the AGN contribution 
to the ionizing background peaks at 2: 2. These results 
are lower than previous estimates of the AGN contribu- 
tion, confirming that ionization from AGNs is insufhcient 
to maintain the observed ionization of the intergalactic 
medium (IGM) at z > 3 (see Bolton et al. 2005 and 
Meiksin 2005 for earlier estimates that suggested this re- 
suh). 

Ideally, however, one would like a direct measurement 
of the shape of the input ionizing source spectrum. For- 
tunately, with the advent of the Galaxy Evolution Ex- 
plorer (GALEX) satellite (Martin et al. 2005), this is 
now possible at 2 ~ 1. Obtaining the shape of this 
spectrum is the primary goal of the present paper. We 
use GALEX grism spectroscopic observations to mea- 
sure the AGN contribution to the ionizing background 
at 2 ^ 1 directly from the UV spectra. For this we use 
the Lya selected sample from Cowie et al. (2010; here- 
after, CBHIO), which should provide an essentially com- 
plete sample of broad- line AGNs (BL AGNs) and a nearly 
complete sample of all AGNs contributing to the ioniz- 
ing flux in the redshift interval 2 = 0.65 — 1.25. We then 
go a step further and measure the ionizing source spec- 
trum at rest-frame wavelengths 700 — 900 A. We note, 
however, that this wavelength coverage is insufficient to 
characterize the full Lyman continuum background. 

The outline of the paper is as follows. In Section [2] we 
briefly describe the 139 2 = 0.65 — 1.25 sources in the 
CBHIO sample. In Section|3]we test the completeness of 
the sample by comparing with an X-ray selected AGN 
sample in the same redshift interval. In Section 2] we 
construct the rest- frame UV (1450 A) luminosity func- 
tion and compare it with the UV luminosity functions 
derived from X-ray and optically selected samples in the 
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same redshift interval. In Section[S]we determine the ion- 
izing source spectrum at z 1 and show that the normal- 
ization agrees well with recent determinations at this red- 
shift. In Section[6]we compute the shape of the comoving 
volume cmissivity (after correcting for incompleteness 
and for the effects of Lyman scattering and absorption in 
the intervening intergalactic medium) versus rest-frame 
wavelength for the redshift intervals z — 0.65 — 0.95 and 
z = 0.95 — 1.25. In Section[7]we summarize our results. 

We use a standard i/o = 70 km s~^ Mpc~^, = 0.3, 
and r^A = 0.7 cosmology throughout. All magnitudes 
are given in the AB magnitude system, where an AB 
magnitude is defined by tuab = — 2. 5 log — 48.60. Here 
is the fiux of the source in units of erg cm"^ s^^ Hz"^ 

2. THE GALEX SAMPLE 

CBHIO constructed a sample of 139 Lya emission line- 
selected sources in the redshift interval z = 0.65 — 1.25 
from nine blank high galactic latitude fields with the 
deepest GALEX grism observations. They took the 
near-UV (NUV) and far-UV (FUV) extracted spectra 
and the NUV magnitudes from the Multimission Archive 
at STScI (MAST). Morrissey et al. (2007) detail the spec- 
tral extraction techniques used by the GALEX team in 
analyzing the grism data and discuss the UV spectra and 
magnitudes. We summarize the fields in Table [1] where 
we give in column (1) the GALEX name, in columns (2) 
and (3) the J2000 right ascension and declination, in col- 
umn (4) the exposure time in kiloseconds, in column (5) 
the limiting NUV magnitude to which the spectra were 
extracted, in column (6) the galactic E{B — V) in the 
direction of the field from Schlegel et al. (1998), in col- 
umn (7) the galactic latitude, in column (8) the number 
of sources with spectra lying within a radius of 32f5 from 
the field center (chosen because in the outermost regions 
of the fields there is a higher fraction of poor quality 
spectra), and in column (9) the number of these sources 
found to have a Lya emission line in the redshift range 
z = 0.65 — 1.25 based on their NUV spectra. The ex- 
tracted sources per field constitute nearly complete sam- 
ples to the NUV limiting magnitudes given in the table 
(i.e., these are NUV selected samples and comprise both 
galaxies and a significant number of stars). The FUV 
spectra cover a wavelength range of - 1300 - 1800 A at 
a resolution of 10 A, and the NUV spectra cover a 
wavelength range of - 1850 - 3000 A at a resolution of 
~ 25 A. 

Each of the spectra was corrected for the Milky Way 
extinction in the field using the Schlegel et al. (1998) 
extinction and the Cardelli et al. (1989) reddening law. 
Because the Galactic extinction towards these fields is 
very low (Table [IJ, the net effect is extremely small. On 
average the shape of the spectra is nearly unchanged and 
the normalization is raised by a factor of 1.08. The low 
extinctions in the present sample represent a consider- 
able advantage over random quasar samples where the 
extinction for many of the objects will be high. 

CBHIO determined the observed area as a function of 
NUV magnitude for each field by computing the ratio 
of sources with GALEX spectra at a given NUV mag- 
nitude to sources with that NUV magnitude in the con- 
tinuum catalog. They multiplied this ratio for each field 
by the area corresponding to the 32. '5 selection radius 



and summed the results for the 9 fields to form an area- 
magnitude relation (their Figure 4(c)). The area is the 
8.2 deg^ area of the 9 fields for magnitudes brighter than 
NUV= 20.5 and then drops as the limiting magnitudes 
of each of the individual fields is reached. It falls to zero 
at ^ 21.8, which is the limiting depth of the deepest 
(GROTH 00) field. 

CBHIO describe their search procedure for finding 
emission lines in the GALEX spectra. Because the spec- 
tra become very noisy at the edges of the wavelength 
ranges, they restricted their higher redshift search to the 
redshift interval z = 0.65 — 1.25 corresponding to the 
wavelength range 2006 - 2736 A in the NUV spectrum for 
the Lya emission line. They measured redshifts for each 
source with an emission line in this wavelength range and 
split the sources into either an AGN class (if there were 
high-excitation lines; usually OVI, NV, or CIV would 
be present in addition to the Lya line) or a candidate 
Galaxy class (if there were only a single line visible; in 
this case the line was assumed to be Lya for determin- 
ing a redshift, which they showed was an extremely re- 
liable assumption). CBHIO found that nearly all of the 
more luminous sources are AGNs (their Figures 4(a) and 
4(b)), leaving relatively few candidate Galaxies in the 
redshift interval z = 0.65 — 1.25 (see also Deharveng et 
al. 2008). Of these, many were detected in the FUV 
(their Figure 4(b)), suggesting that they are AGNs, too, 
since z ~ 1 galaxies should be very faint or undetected in 
the FUV band, which lies below the Lyman continuum 
edge at these redshifts (e.g., Siana et al. 2007; CBT09). 
We shall hereafter refer to all of the sources that were not 
classified as AGNs as unclassified sources. Since nearly 
all of them are likely AGNs, we will refer to the com- 
bined AGN plus unclassified source sample as our "full" 
AGN sample. We note that we have confirmed four of 
the unclassified sources as star formers through optical 
spectroscopic observations, but excluding them from our 
analysis makes no difference to the results. 

3. TESTING THE COMPLETENESS OF THE 
GALEX AGN SAMPLE 

We test the completeness of the GALEX AGN sample 
by comparing with an X-ray selected AGN sample in the 
same redshift interval. In Figure [T] we show the number 
counts for the GALEX AGN sample (black squares) and 
for the GALEX full AGN sample (blue triangles) versus 
NUV magnitude in the redshift interval z = 0.65 — 1.25. 
(There is only a very small difference between the AGN 
and the full AGN results.) In each magnitude bin the 
counts are the sum of the inverse areas for all the sources 
in that bin. The error bars are ± la based on the Pois- 
son errors corresponding to the number of sources in each 
bin. The number counts per square degree per unit mag- 
nitude, N, are well described by the form (blue dashed 
line) 

log N = 0.40 (NUV - 18.4) (1) 

to the limiting NUV magnitude of the grism data of 21.8. 
This functional form implies an equal amount of light in 
each magnitude interval. 

For the X-ray selected AGN sample in the same red- 
shift interval we use the extremely spectroscopically com- 
plete Ghandra samples of Trouille et al. (2008, 2009) and 
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TABLE 1 
GALEX FIELDS 



Name 


R.A. 


Dcd. 


Exp. 


NUV„„ 


E{B - V) 


Gal. Lat. 


No. w/sp. 
(< 32.'5) 


No. w/Lyo 




(J2000.0) 


(J2000.0) 


(ks) 


(mag) 




(deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


GROTH 00 


214.99182 


52.78173 


291 


21.8 


0.007 


59.5 


1152 


61 


NGPDWS 00 


219.15610 


35.17135 


156 


21.5 


0.009 


66.2 


738 


39 


CDFS 00 


53.12779 


-27.87137 


149 


21.5 


0.008 


-54.4 


876 


36 


COSMOS 00 


150.11900 


2.20583 


140 


21.5 


0.018 


42.1 


779 


27 


ELAISSl 00 


9.63857 


-43.99023 


84 


21.2 


0.007 


-72.9 


589 


17 


SIRTFFL 00 


259.12387 


59.90915 


80 


21.2 


0.020 


35.0 


800 


21 


LOCK 00 


162.67843 


58.73117 


48 


20.9 


0.009 


52.1 


390 


10 


SIRTFFL 01 


260.41425 


59.34286 


34 


20.7 


0.029 


34.4 


612 


14 


HDFN 00 


189.20946 


62.19772 


24 


20.5 


0.012 


54.8 


265 


10 
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Fig. 1. — Number counts of the GALEX AGNs in the redshift 
interval z = 0.65 — 1.25 vs. NUV magnitude (black squares with 
1(T errors based on the number of sources in each bin). The blue 
triangles show the values if we also include the unclassified sources 
in this redshift interval. For comparison, we show the number 
counts of an X-ray selected sample in a 0.9 deg^ area (the CLANS 
and CLASXS fields) with red diamonds and in a smaller but much 
deeper 200 arcmin^ area (the CDF-N) with purple diamonds. Over 
the magnitude range of the GALEX grism data the GALEX AGN 
number counts agree well with the X-ray selected AGN number 
counts and are well fit with a slope of 0.4 (blue dashed line). Thus, 
each magnitude bin contains roughly equal light density. 



references therein on the Chandra Large- Area Synopotic 
X-ray Survey (CLASXS), the Chandra Lockman Area 
North Survey (CLANS), and the Chandra Deep Field- 
North (CDF-N). Essentially all of the X-ray sources in 
the redshift interval z — 0.65 — 1.25 are expected to be 
identified. 

The CLASXS and CLANS fields are covered by a num- 
ber of GALEX pointings, all of which are roughly 30 ks 
in depth. The CDF-N GALEX exposures are 165 ks in 
the NUV and 96 ks in the FUV. CBT09 describe how the 
GALEX magnitudes for the sources in the X-ray sam- 
ples were measured. In brief, given the large point-spread 
function (PSF) of GALEX (4'.'5-6'.'0 FWHM), they used 
an 8" diameter aperture to measure the magnitudes us- 
ing the GALEX zeropoints from Morrissey et al. (2007), 
and they corrected these to approximate total magni- 
tudes using an offset that they determined to be —0.41. 
They measured the noise level in the images by randomly 
positioning apertures on blank regions of sky and mea- 



suring the dispersion. In the CLASXS and CLANS fields 
they found the la limits to be slightly variable but typ- 
ically about 26.1 in the NUV images and 26.5 in the 
FUV images. In the CDF-N field they found la limits 
of 26.8 in the NUV image and 27.4 in the FUV image. 
Only a small fraction (^ 4%) of the X-ray sources in the 
CLASXS or CLANS fields are not covered by one of the 
GALEX pointings, and all of the CDF-N sources are cov- 
ered. CBT09 only eliminated four sources where visual 
inspection showed that there was contamination from a 
nearby brighter GALEX source at the same wavelength. 
The X-ray sources themselves are so sparse that contam- 
ination by another X-ray source may be neglected. 

In Figure [1] we show separately the number counts 
from the shallower CLASXS and CLANS fields (red 
diamonds), where the 2 — 8 keV flux limit is ~ 7 x 
2Q-15 gj-g cjj^-s g~i g^j^j^ ^Yie area is ~ 0.9 deg^, and 
the number counts from the CDF-N (purple diamonds), 
which is about 50 times deeper but has a much smaller 
area of ^ 200 arcmin^. Over the magnitude range 
of the GALEX grism data the GALEX AGN num- 
ber counts agree extremely well with the CLASXS and 
CLANS AGN number counts. When we compare the 
NUV counts from the shallower X-ray fields with those 
from the CDF-N, we see that the former are becoming in- 
complete fainter than NUV ~ 22.5. The CDF-N number 
counts themselves should be complete to several magni- 
tudes fainter than this, so their behavior suggests that 
the AGN number counts begin to turn down from the 
0.4 slope at NUV ^ 23.5. This in turn suggests that the 
GALEX AGN sample measures a large fraction of the 
total NUV light at these redshifts. 

4. THE AGN LUMINOSITY FUNCTION 

We next computed the binned rest-frame UV (1450 A) 
luminosity function for the GALEX full AGN sam- 
ple using the NUV-continuum magnitudes and the 1/V 
method (Felten 1976). We made the small differential 
i^T- corrections to 1450 A by assuming a quasar spec- 



tral energy distribution (SED) with 



(Vanden 



Berk et al. 2001). We show this luminosity function with 
black squares in Figure [21 The error bars are ±lcr based 
on the Poisson errors corresponding to the number of 
sources in each bin. We compare with the 1450 A lumi- 
nosity function derived from the X-ray selected BLAGN 
luminosity function of Yencho et al. (2009) (red dia- 
monds). For the conversion of 2 — 8 keV luminosity 
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Fig. 2. — The rest-frame UV (1450 A) luminosity function of the 
GALEX full AGN sample in the redshift interval z = 0.65 - 1.25 
(black squares with la error bars from the observed number of 
sources in each bin). Red diamonds show the converted X-ray 
selected BLAGN luminosity function in the same redshift in- 
terval from Yencho et al. (2009). The green curves (lower for 
2 = 0.68 — 0.95; upper for z = 0.95 — 1.25) show the converted 
optically selected luminosity functions from Groom et al. (2004). 

(Lx) to 1450 A absolute magnitude we used the rela- 
tion iyL^{2300 k)/Lx = 7.46 x {Lx/IG^^ erg s'^f-^^ 
from CBT09, which is similar to previous work (e.g., Vi- 
gnali et al. 2003). We again made the small differential 
corrections to 1450 A by assuming ~ v^^-^^. Fi- 
nally, we compare with the optically selected quasar lu- 
minosity function of Croom et al. (2004) (green curves) 
scaled to 1450 A by assuming f^, ^ v^^-^^. A similar 
result is obtained using the Richards et al. (2005) lumi- 
nosity function. 

The GALEX luminosity function agrees extremely well 
with these other determinations and reaches about a 
magnitude fainter than the wide-field optical samples 
used in the Croom et al. (2004) measurement. The X- 
ray-determined luminosity function from Yencho et al. 
(2009) is deeper yet and shows a turn-down at fainter 
magnitudes. This means that the UV light density is 
nearly convergent in the GALEX full AGN sample at 
these redshifts. If we extrapolate the GALEX luminos- 
ity function with a constant value below A/(1450)ab = 
—22, then the additional UV light density at magnitudes 
fainter than M(1450)ab = -22 is only 24% of the total 
light for the observed sources with magnitudes brighter 
than this value. 

5. THE IONIZING SPECTRUM OF Z - 1 AGNS 

We next compared the flux measured from the spec- 
tra at rest-frame wavelengths both above and below the 
Lyman continuum break with the selection NUV mag- 
nitude. In Figure |3] we show the observed flux from 
the GALEX spectra at rest-frame (a) 800 - 880 A (for 
z = 0.65 — 1.15 only; note that we use only the portion 
of the wavelength range that is covered by the GALEX 
FUV spectra) and (b) 1400 - 1500 A (for z = 0.65 - 0.95 
only, in order for the GALEX NUV spectra to cover 
the rest-frame wavelengths) versus NUV. The tight re- 
lation in (b) simply reflects the fact that the redshifted 
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Fig. 3. — (a) Observed flux from the GALEX spectra at rest- 
frame 800 — 880 A vs. NUV magnitude for sources in the redshift 
interval z = 0.65 — 1.15. Red squares denote sources classified 
as AGNs due to the presence of high-excitation lines, while black 
squares show the unclassified sources in this redshift interval. The 
blue line shows a linear fit to the data, (b) Same as (a) but for 
rest-frame 1400 — 1500 A and restricted to z = 0.65 — 0.95, where 
the GALEX NUV spectra cover the rest-frame wavelengths. The 
blue line shows the expected value of the flux determined from 
the NUV continuum magnitude, since redshifted 1400 — 1500 A 
matches closely to the NUV wavelength. There is a very small 
multiplicative factor (0.95) between the flux measured from the 
spectra and the flux measured from the broadband magnitudes. 

1400 — 1500 A spectral interval matches closely to the 
NUV wavelength. The data provide a test of the GALEX 
spectral intensity calibration, which appears to be ex- 
tremely accurate, since the spectral fluxes differ by less 
than 5% from the broadband continuum fluxes. 

As can be seen from Figure Ela) , most of the sources 
also obey a linear relation between the spectral flux at 
wavelengths just below the continuum break and the 
NUV flux. At 840 A we find 



/a(840 A) = 1.65 X 10-16 io-o.4(Nuv-20) 



(2) 



where f\{8A0 A) is the average spectral flux of the sources 
at rest-frame 800 — 880 A in erg cm~^ A"^ . We show 
this as the blue line in Figure[3Ua). A small fraction of the 
sources (i.e., those where there is a signiflcant continuum 
break) lie significantly below this relation. The mean fx 
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Fig. 4. — Cumulative comoving volume emissivity (VEM) for 
the full AGN sample at rest-frame 840 A vs. rest-frame 1450 A 
absolute magnitude for the two redshift intervals 2 = 0.65 — 0.95 
(red curve) and z = 0.95 — 1.15 (black curve). The horizontal 
lines show the incompleteness corrected values for the two redshift 
intervals. 

ratio, /a(840 A)//a(1450 A), is 1.11; the mean ratio is 
0.37. CBT09 gave an U ratio of 0.26 between 912 A and 
2300 A, which for - ly-^-^^ (Vanden Berk et al. 2001) 
would correspond to 0.34 for the f^{840 A)//^(1450 A) 
ratio, in close agreement with the present value. 

We then computed the comoving volume emissivity 
(VEM) in the rest-frame wavelength range 800 - 880 A 
using the 1/V methodology. In Figure U] we show the 
cumulative contribution to the VEM as a function of the 
rest-frame 1450 A absolute magnitude for the redshift in- 
tervals z = 0.65 — 0.95 and z = 0.95 — 1.15, which reach 
Mab(1450 A) 21 and 21.7, respectively. 

We have taken a simple approach to estimate the in- 
completeness correction due to sources fainter than our 
selection NUV magnitude limit. As a result of the two 
linear relations in Figure O the /a (840 A)//a(1450 A) 
ratio stays fixed over the NUV magnitude range of our 
sources, which implies that the spectral shape is invariant 
over that magnitude range. If we assume that the ratio 
also remains fixed at fainter magnitudes, then any esti- 
mate of the incompleteness correction made at 1450 A 
can be applied at any other wavelength. 

In order to estimate the incompleteness correction at 
1450 A, we used the form of the 1450 A luminosity func- 
tion (e.g., Figure[5]). We determined the amount of miss- 
ing hght at 1450 A by assuming a constant value per 
Mpc"^ per mag below the limiting absolute magnitude 
that we measured. From the ratio of this integral to the 
integral made only to the limiting absolute magnitude, 
we obtained correction factors of 1.09 (z — 0.65 — 0.95) 
and 1.24 (z — 0.95 — 1.15). We applied these corrections 
to the measured 840 A VEMs, giving incompleteness cor- 
rected values that are not substantially higher than the 
direct measurements (horizontal lines in Figure |4]). Nu- 
merically these values are 1.07±0.19x 10^'* erg Mpc~^ s~^ 
at z = 0.8 and 1.68 ± 0.23 x lO^^ erg Mpc'^ at 
z = 1.05. 

We can compare these results to CBT09's determina- 
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Fig. 5. — Evolution of the comoving volume emissivity at rest- 
frame 840 A vs. redshift for the redshift intervals z = 0.65 — 0.95 
and 2 = 0.95 — 1.15 (blue open squares show the directly mea- 
sured values, and solid squares show the incompleteness corrected 
values). We also show the evolution determined by CBT09 from 
measuring the GALEX FUV magnitudes of an X-ray selected sam- 
ple (red diamonds). For each point the x-axis error bars show the 
redshift interval and the y-axis error bars show the Icr error com- 
puted using the jackknifc method. 

tion of the redshift evolution of the VEM just below the 
Lyman continuum edge. In Figure [5] we show CBT09's 
X-ray based points with red diamonds and the present 
points with blue squares (open for directly measured 
values; solid for incompleteness corrected values). The 
Icr errors for both samples have been computed with 
the jackknifc method. The smaller error bars on the 
present points reflect the much larger area sampled by 
the GALEX images than by the X-ray images. It is very 
reassuring that the two samples agree within the formal 
errors. 

6. THE SHAPE OF THE SPECTRUM 

With the grism data we are able to move one criti- 
cal step beyond CBT09 and compute the shape of the 
VEM with wavelength. As discussed in Section [1] a 
direct measurement of the shape of the input ionizing 
source spectrum is ideally required for the cosmological 
modeling of galaxies and the intergalactic gas. For each 
rest-frame wavelength interval we computed the VEM 
using the same procedure that we used for the rest-frame 
800 — 880 A interval, including applying at all wave- 
lengths the incompleteness corrections of 1.09 and 1.24, 
respectively, for the low and high-redshift intervals. How- 
ever, in order to obtain the true rest-frame VEM spec- 
trum, we also need to correct for the effects of the Lyman 
scattering and absorption in the intervening intergalac- 
tic medium (M0ller & Jakobsen 1990; Haardt & Madau 
1996; Fardal et al. 1998). While this correction is not 
large at z = 1, it does slightly steepen the spectral slope 
in the extreme UV (EUV) and raise the normalization. 

The effects of the intergalactic medium (IGM) may 
be separated into scattering by the Lya forest lines be- 
low a rest-frame of 1216 A and photoelectric absorption 
from the Lyman continuum of hydrogen below 912 A. 
At the wavelengths of interest we may neglect helium. 
The effects of the forest are dominated by systems with 
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column densities around 10"'^^ cm^'^, where Lya lines sat- 
urate, while the photoelectric absorption is dominated by 
systems with column densities ~ 10^^ cm~^, where the 
Lyman continuum edge saturates. 

The decrement in the Lya-only portion of the forest 
has been directly measured by Kirkman et al. (2007). 
They obtained a value of 0.037 x (1 + z)^-^ at these red- 
shifts. The redshift dependence, which causes less light 
to be absorbed at shorter wavelengths (i.e., because there 
are fewer systems at lower redshifts) , is roughly balanced 
by the onset of absorption by the higher Lyman series 
lines at shorter wavelengths, resulting in a roughly con- 
stant reduction as a function of wavelength near z 1 
(see Figure 6a of Inoue & Iwata 2008). 

The photoelectric absorption is most directly related 
to the density of Lyman limit systems (LLSs). In LLSs 
the optical depth, t, at the Lyman continuum edge ex- 
ceeds one. The surface density of LLSs at these red- 
shifts has recently been remeasured by Songaila & Cowie 
(2010). They found that the number density per unit 
redshift evolves as 0.15 x (1 + z)^'^. This normalization 
is slightly lower than previous measurements at these 
redshifts (Stengler-Larrea et al. 1995). We may use this 
normalization in conjunction with the distribution of the 
number of systems as a function of the neutral hydrogen 
column density to compute the complete photoelectric 
effect. The column density distribution is usually pa- 
rameterized as a power law or broken power law function 
of the column density. Inoue & Iwata (2008) assumed a 
broken power law with a slope of —1.6 below the Ly- 
man limit and —1.3 above (e.g., Janknecht et al. 2006; 
Songaila & Cowie 2010), which we adopt here. We have 
also computed the correction for a single power law with 
a slope of —1.4 in the vicinity of the Lyman limit (Misawa 
et al. 2007) to investigate the sensitivity to the assumed 
column density distribution. This slightly decreases the 
correction, but the changes produced are smaller than 
the statistical errors for both the normalization and the 
spectral slope. 

We have used these values to make Monte Carlo es- 
timates of the photoelectric absorption (Bershady et al. 
1999; Inoue & Iwata 2008). In Figure [B] we show the 
photoelectric absorption decrements that we obtained at 
(a) z = 0.8 and (b) z = 1.1 from averaging (in each case) 
1000 lines of sight. We show the values obtained both 
from using only systems with optical depths less than 
one at the Lyman edge (dashed line) and from using all 
the systems regardless of their optical depths (solid line). 
The Monte Carlo simulations also allow us to estimate 
the uncertainty introduced by the finite number of sys- 
tems we are observing. For the case where we used all the 
systems regardless of their optical depth, we also show 
with dotted lines the range of corrections that we found 
in 20 computations using only 60 lines of sight (which 
is roughly equal to the number of lines of sight in the 
present sample; we have 53 systems at z = 0.65 — 0.95 
and 63 at z = 0.95 — 1.25). At the shorter wavelengths 
this corresponds to about a 10% uncertainty in the cor- 
rection. The large range is a consequence of the small 
number of systems which contribute to the opacity at 
T > 1. The uncertainty in the correction for the r < 1 
systems only is much smaller (about 3% or less). This 
is because there is a much larger number of low column 



density systems, and hence the correction is better aver- 
aged even in a small number of lines of sight. 

We have applied the correction in two ways. In the first 
we excluded the small number of quasars with measured 
LLSs in their spectra using the list in Songaila & Cowie 
(2010). We corrected for the missing light from these 
AGNs by proportionally reducing the observed area. We 
then used the r < 1 correction to compute the final VEM. 
In Figure [iKa) we show the VEM corrected in this way 
and its la errors in the redshift intervals z = 0.65 — 0.95 
(red squares) and z = 0.95 — 1.25 (blue diamonds) versus 
wavelength. The advantage of this method is that it 
minimizes the correction and also avoids the statistical 
uncertainties in the r > 1 correction. However, we may 
fail to identify LLSs in the fainter and noisier spectra. 
We have therefore also computed the VEM from the full 
quasar sample (i.e., not excluding systems with LLSs). 
In this case we apply the full correction, including the 
opacity from the r > 1 systems. In Figure [Tl^b) we show 
the VEM corrected in this way and its la errors in the 
same redshift intervals. We summarize the values that 
we derived using both methods in Table 2. 

Scott et al. (2004) fitted a single power law to their 
composite EUV QSO spectrum from Far- Ultraviolet 
Spectroscopic Explorer (FUSE) data, but both Zheng et 
al. (1997) and Teller et al. (2002) used a broken power 
law to fit their composite EUV QSO spectra from Hubble 
Space Telescope (HST) data. For each redshift interval 
we have fitted a broken power law to the data, exclud- 
ing the bins at 775 A, 1025 A, and 1225 A to exclude 
the strong emission lines in the spectra. The position 
of the break is quite uncertain. We have chosen to po- 
sition it at 1225 A. We show the fits as straight lines 
in Figure [7] for each case. Although it is possible that 
there is more structure within the error bars, the present 
data do not justify a more complicated fit. We present 
this simple parameterization in order to provide quanti- 
tative numbers for use by others. Since the differences 
between the two methods lie within the statistical errors 
(see Table 2), we consider only the cases with the LLSs 
excluded (Figure [7^), which should be more robust. For 
z = 0.65 — 0.95 we obtain 

log(VEM) = (24.11 ± 0.03) -I- (1.76 ± 0.25) log(w/800 A) 

(3) 

below 1225 A and 

log(VEM) = (24.46±0.01) + (0.36±0.45)log(?i;/1250 A) 

(4) 

above 1225 A. For z = 0.95 - 1.25 we obtain 

log(VEM) = (24.34 ± 0.02) + (2.01 ± 0.20) log(w/800 A) 

(5) 

below 1225 A and 

log(VEM) = (24.69±0.02)-(0.51±0.97)log(w/1250 A) 

(6) 

above 1225 A. In all of the above equations, w is the 
rest-frame wavelength. 

The slopes are consistent between the two redshift 
ranges within the statistical errors, while the normal- 
ization rises by 0.23 dex or a factor of 1.7. The short 
wavelength slope is considerably steeper than the value 
of 0.56 (-0.28, 0.38) derived by Scott et al. (2004) but 
agrees well with the values of 1.76 ±0.1 2 derived by Teller 
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Fig. 6. — Monte Carlo estimates of the photoelectric absorption 
decrements at (a) z = 0.8 and (b) z = 1.1 obtained from averaging 
(in each case) 1000 lines of sight. The dashed curves only show 
the values from systems with optical depths less than one at the 
Lyman edge, while the solid curves show the values from all systems 
regardless of their optical depths. The dotted curves show the 
range of corrections for 20 computations of only 60 lines of sight, 
which is similar to the number of systems that we have in each 
redshift interval. 

et al. (2002) and 1.96 ± 0.15 derived by Zheng et aL 
(1997). These various samples have different himinos- 
ity and redshift selections, which may account for some 
part of the differences. The present sample is properly 
weighted over the sources contributing to the ionizing 
background. 

7. SUMMARY 

We used a complete sample of Lya emission-line se- 
lected AGNs obtained from nine deep blank fields ob- 
served with the grism spectrographs on GALEX to mea- 
sure the AGN contribution to the ionizing background 
at z 1. We found that the GALEX AGN luminos- 
ity function is in good agreement with luminosity func- 
tions obtained from optical and X-ray AGN samples. We 
also found that our measurements of the 840 A VEM at 
z = 0.8 and z = 1.05 agree within the formal errors with 
the redshift evolution of the comoving volume emissivity 
determined by CBT09 using an X-ray selected AGN sam- 
ple. Finally, for the first time, we were able to compute 



Fig. 7. — Shape of the comoving volume emissivity vs. rest-frame 
wavelength for the redshift intervals z = 0.65 — 0.95 (red squares) 
and z = 0.95 — 1.25 (blue diamonds). The values shown are cor- 
rected for incompleteness, and the error bars are ±1 cr computed 
using the jackknife method. In (a) we exclude the quasars with 
measured LLSs in their spectra before applying the t < 1 cor- 
rection to compute the final VEM. In (b) we use the full quasar 
sample and apply the full correction, including the opacity from 
the T > 1 systems, to compute the final VEM. The straight lines 
show broken power law fits. 



directly the shape of the VEM with wavelength for the 
redshift intervals z = 0.65 — 0.95 and z = 0.95 — 1.25. We 
corrected both for incompleteness and for the effects of 
the Lyman scattering and absorption in the intervening 
intergalactic medium. The slopes are consistent between 
the two redshift ranges within the statistical errors, while 
the normalization rises by a factor of 1.7. At < 1125 A 
our slope is considerably steeper than that determined 
by Scott et al. (2004) using FUSE data, but it is in good 
agreement with the slopes determined by Teller et al. 
(2002) and Zheng et al. (1997) using HST data. 
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AST-0708793 (A. J. B.) and AST-0709356 (L. L. C.), 
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* Raw values do not have the incompleteness correction or the Lyman correction. The VEM and errors are given in units of 
10^4 erg Mpc-3 s"! Hz"!. 



